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Abstract 

The first three-dimensional simulation of shear-induced phase transitions in 
a polymeric system has been performed. The method is based on dynamic 
density-functional theory. The pathways between a bicontinuous phase with 
developing gyroid mesostructure and a lamellar/cylinder phase coexistence 
are investigated for a mixture of flexible triblock ABA-copolymer and solvent 
under simple steady shear. 
Various self assembly systems such as lyotropic liquid crystals, surfactants, block copoly- 
mers can form ordered mesophases (lamellar, cylindrical, spherical, etc). These phases have 
received much attention because of the fundamental interest to establish universal laws for 
self-organization phenomena, and also because of the wide range of applications in materials 
scienceS^. 

An interesting issue in the design of new materials is modulation of phase behaviour 
by external and internal factors such as fiowJ^, reactionsi, temperature inhomogeneityi, 
confinements and surfaces^. In particular, externally applied shear flows are found to 
lead to macroscale order in block copolymer systems. Moreover, shear introduces a new 
kind of phase behaviour of block copolymer systems, the so-called orientational phase 



transit ions'^^^^ 



So far, lamellar, hexagonal cylindrical and cubic micellar phases of block copolymers 



under shear have been thoroughly investigated in experimentlli§'[i'El. Theoretical under- 
standing is partly reached for the lamellar geometry0lil"El and for the hexagonal cylindrical 
phase0'Ei0. Experimental observation of a sponge phase in low molecular weight surfac- 
tants under shear has been carried out recently ini. Elongation of mesostructures has been 
recently found in a sheared isotropic bicontinuous polymer phased. More complex phases 
such as gyroid and coexistence of phases in shear flow still require experimental examination. 

Previously, computer simulations of polymer morphologies under shear have been carried 
out for lamellar and cylindrical phases in 2D systems (see references in0). Recently we have 
reported 3D density functional calculations for these two phasesiiS Here we report on the 
effect of shear on more complex phases — bicontinuous gyroid and lamellar / cylinder coexis- 
tence in a three component copolymer system. We have observed shear-induced transitions 
between these states. Many shear-induced transitions are found in experiments and theory^. 
Transitions we present here require future experimental and theoretical investigation. 

The time evolution of the density field under simple steady shear flow, Vx = jy, 
Vy = Vz = 0, can be described by a time dependent Landau- Ginzburg type equation with 
a convective ternflilii and a stochastic termilil. In general, the velocity field can be 
found from hydrodynamics. For a system with different viscosities of the components, this 
leads to an adaption of the linear velocity profiled. Such refinement is important for the 
description of temperature dependency of the orientation of mesostructure lattice in the 
gradient- vorticity plane at high shear (see discussion in0). However, in the weak segrega- 
tion regime, when there are no steep concentration gradients, the basic features of the process 
of alignment in flow (regardless of fine details of orientational transitions) can be described 
accurately by diffusion-convection equation, with an imposed velocity profilelllilB. 

In contrast to traditional schemes of polymer phase separation dynamics where a Landau 
Hamiltonian is used with vertex functions calculated following the Random Phase Approx- 
imation (see e.(7.0), we numerically calculate the "exact" free energy F of polymer system 
consisting of Gaussian chains in mean field environment using path integral formalism@0. 
The benefit of our approach is that it avoids the truncation of the free energy and therefore 



represents intermediate (metastable) states more accurately. 

Some time ago, Matsen and Schickl introduced a powerful method for SCF calculations 
of equilibrium block copolymer morphologies. Our approach uses essentially the same free 
energy functional and complements the static SCF calculations by providing a dynamical 
picture of the system which is crucial for systems under sheari'SS. In our method the 
calculations are in direct space, without any bias with respect to the morphology formation. 

The system presented here is a mixture of Gaussian chains E^PgE^ (bead names are 
arbitrary) and solvent which is parameterized to model 60% aqueous solution of a triblock 
copolymer Pluronic surfactant PL64 {EO)i3{PO)3oiEO)i3 with the hydrophobic block in 
the mi ddle§. The choice is justified by huge variety of available experimental data. In 
experimental phase diagram, 60% polymer concentration corresponds to a very complex 
phase coexistence region^. 

Fig. 1 gives a schematic overview of the simulation of the system in a 3D box (64x64x64). 
The starting configuration is a homogeneous distribution of the components. The system 
(Fig. 0) demonstrates the development of a bicontinuous morphology with clear gyroid-type 
connectivity (Fig. however still without global symmetry throughout the sample. The 
system can remain in this phase very long, slowly rearranging the structure and keeping 
gyroid-type connectivity (see also Fig. ^). 

After applying shear (7 = 10~^) to the morphology shown in Fig. |^, the system slowly 
deforms. However, the connectivity hardly changes: the system remains bicontinuous and 
similar to the morphology shown in Fig. ^ even at large shear strains 7 > 1(100%). While 
shearing continues, connections in the system start to break, separate pieces of the structure 
reconnect to form a pattern aligned in the flow direction. New structure consists of coexisting 
lamellae and hexagonally packed cylindrical clusters, as shown in Figs, ^a and |a. 

Stopping the shear at some point (r = 7500 (7 = 11), Fig. |1|) leads to the reorganization 
of the structure via migration of defects (holes in lamellae, necks between cylinders). Fig. ^. 
The time evolution of deformed and broken bicontinuous structure involves a lot of perforated 
lamellar clusters. Migration of defects in these clusters leads sometimes to intermediate 
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structures with locally hexagonal-like arrangement of holes (Fig. middle picture). After 
long relaxation without shear the structure remains a coexistence of defected cylinders and 
lamellae (Fig. ^) but with a better ordering compared to the moment of stopping shear 
(Fig. ^a). The free energy plot (Fig. |^) provides information which state is metastable. 
Bicontinuous and coexistence states are very close in free energy value in plateau region 
(Fig. but the last is slightly lower. 

A second period of stronger shear (7 = 5- 10~^) very fast breaks up the remaining con- 
nections and the system flows as a lamellar /cylinder coexistence without changing lamel- 
lar/cylinder volume ratio while shear continues, r = 23500 (7 = 35), Fig. The lamellae 
have only a few holes. The cylinder region no longer has "neck" defects and the only re- 
maining defect is a long living dislocation which remains stable during the whole period of 
shearing. As a result the hexagonally packed cylinders consist of two clusters with different 
internal hexagonal orientation (as can also be seen from the double peaks in the structure 
factor. Fig. ^). After stopping shear at 7 = 35 (Fig. ^) the system relaxes to a lamel- 
lar/cylinder coexistence in which the hexagonal lattice has only one orientation (Figs, ^b, 
^ and 1^). Undulated cylinders partly transform into lamellae via formation of subsequent 
necks. Fig. |^. The lamellar phase still consists of perforated lamellae, with a low fraction of 
holes. This coexistence of phases seems to be stable and even application of a much higher 
noise {Q = 100, the region is schematically shown in Fig. [^) does not change the picture 
considerably. 

Both periods of shear demonstrate that the lamellar cluster becomes larger after shear is 
released (Figs. |^,b and |c,d). This is consistent with the fact that switching on the shear 
again squeezes the size of the lamellar region. Figs. ^ and 

The stability of the phase coexistence morphology. Figs. is challenged by applying 

shear (7 = 5- 10~^) in the z-direction (perpendicular to lamellae/cylinders). Fig. |1]. In Figs. 



TO| a and |TO|c we see the result of shearing for a period of r = 50 (7 = 0.5) and r = 300 
(7 = 3), respectively. While shearing the lamellae start to tilt and thin. The cylinders start 
to roll over each other and deform to a slightly prolate shape in cross-section. This leads to 



an energetically unfavorable cubic cylindrical lattice, Fig. |T^a, or to a morphology with very 
oblong cylinders in cross-section and partly broken lamellae, Figs. p!0|c,p]l|a. Releasing shear 
at these different stages leads to different phenomena. After a relatively small distortion 
of the system (7 = 0.5, Fig. p!0|a) in short time the system relaxes back from cubic to 
hexagonal packed cylinders, keeping the lamellar cluster intact. Fig. pUjb. However, if shear is 
stopped after a larger distortion (7 = 3, Fig. [lO|c), then very fast necks form throughout the 
sample, Fig.[To|d. The structure becomes bicontinuous with sometimes obviously gyroid-like 
connectivity. Fig. |Tl| c. However the global structure is different from the initial morphology 
in Fig.|^. The process of relaxation after stopping shear goes via following stages: frustrated 
lamellae/cylinders bicontinuous phase lamellae/cylinders. The first stage is relatively 
fast whereas the last is slow, which can be viewed in hardly noticeable changes in structure in 
Figs.pJ]c and d. Free energy plot (Fig. ^) also demonstrates that system tends to relax back 
to a coexistence state. Thus, in this case the bicontinuous phase is a long living intermediate 
stage with the free energy very close to lamellae/cylinders coexistence state. 

In summary, we have performed the first 3D shear simulation of a bicontinuous ABA- 
copolymer /solvent system. New shear induced phase transitions from bicontinuous phase 
to lamellae/cylinder coexistence and back have been detected. These two states have very 
close free energy values and can be separated by a barrier in the free energy landscape. As 
a result the polymeric system can be trapped in either of them. 

A.V.Z. and G.J.A.S. acknowledge support of the MesoDyn project ESPRIT No. EP22685 
of the European Community. We acknowledge support of NCF (Stichting Nationale Com- 
puter Faciliteiten) . 
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FIG. 1. The sample history. Arrows indicate figures with snapshots. The box with the curve 
indicates a high noise region and the symbol _L indicates change in shear direction. 7 = At'j is 
the dimensionless shear rate and r is dimensionless time with time step At = (3~^Mh^'^At = 0.5 
(seei). 




FIG. 2. An isosurface of PO blocks at r = 2000 (no shear). 




FIG. 3. Evolution of a detail of the bicontinuous morphology in the absence of shear at 
r = 2000, 4975, 7500 (from left to right). 
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FIG. 4. Morphology of PO blocks at r = 7500 (a) and 50000 (b). 




FIG. 5. Orthogonal projection in x-direction of the morphologies at different times: r = 7500 
(a), 20000 (b), 23500 (c), 50000 (d). The arrow in (c) denotes a dislocation. 
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FIG. 6. Detail of a perforated lamellae transformation: r=7500, 9500, 18500 (from left to right). 
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FIG. 7. Free energy as function of the time for the first (a) and the last (b) period of shear. 
The inset is magnified representation of the left plot. Arrows indicate figures with snapshots. First 
and second arrows from the left correspond to starting end stopping times of shearing (c/. Fig]l|). 
The dashed line on the left graph corresponds to bicontinuous phase evolution in the absence of 
shear (see Fig.^). 




FIG. 9. 3D structure factor summed in x-direction (log-scale) for r = 23500 (a), (c/. Fig.5c), 
and r = 50000 (b), (c/. Fig.5d). 
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FIG. 10. Results of shear in z-direction, taking the morphology as shown in Fig.^b ( r = 50000), 
as starting structure. Orthogonal projections of the morphologies at r = 50050 (a), the end of 
shear, and r = 54050 (b). A longer period of shear from the same starting structure: r = 50300 
(c), the end of shear, and at r = 52300 (d). 




FIG. 11. Detail of structure relaxation after shearing in z-direction at r = 50300 (a) (c/. 
Figjgc), 52300 (b) (c/. Fig.^d), 54000 (c), 56500 (d). 
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